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A new type of external fields violating the particle number preservation is studied in one- 
dimensional strongly correlated systems by the Density Matrix Renormalization Group method. 
Due to the U(l) symmetry breaking, the ground state has fluctuation of the total particle num- 
ber, which implies injection of electrons and holes from out of the chain. This charge fluctuation 
can be relevant even at half-filling because the particle-hole symmetry is preserved with the finite 
effective field. In addition, we discuss a quantum phase transition obtained by considering the 
symmetry-breaking fields as a mean field of interchain-hopping. 
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I. INTRODUCTION 

Doped Mott insulators is one of possible candidates of 
superconductors with electron-electron correlation. At a 
rational filling, strong electron-electron interaction makes 
electrons localized in real space. This is the Mott insu- 
lator where charge excitations are gapped. Even at this 
rational filling, the spin degree of freedom survives as a 
gapless mode where quantum objects as collection of the 
S = 1/2 spins form a singlet ground state. 

When mobile carriers are introduced into the Mott in- 
sulator, we may expect that the charge gap is destroyed 
which realizes the superconducting ground state driven 
by the electron interaction. In the resonating valence 
bond (RVB) picture proposed by Anderson to describe 
high-T c superconducting cuplatesQ,Q] , doped holes itin- 
erate in a spin-singlet ground state and condense into a 
superconducting state. If there is no doped hole, i.e., 
half-filled case, the spin-singlet ground state is expected 
as a Mott insulator. 

Apart from real doping, that is, changing chemical po- 
tential, there can be several possibilities for the effective 
carrier doping. One of a natural possibility can be geo- 
metrical frustration in layered organic superconductors, 
which is approximated by a half-filled Hubbard model 
with next-nearest-neighbor hopping|3|. Another one is 
gossamer superconductivity by Laughlin [4( . Even at the 
half-filed case, finite double occupancy may destroy the 
Mott insulator at a small on-site Coulomb repulsion and 
lead theground state to the gossamer superconducting 
state^, H- This theory has been also applied to organic 
superconductors 0- The gossamer superconducting state 
is realized not as a ground state of strongly correlated 
system but a variational function based on a Bardeen 
Cooper Schrieffer (BCS) superconducting state. How- 
ever, it can be an exact ground state of the model Hamil- 
tonian Q which violates the charge conservation as the 
BCS Hamiltonian. It is thus theoretically interesting to 
consider U(l) symmetry breaking generically. 

Let us recall the BCS Hamiltonian. The BCS Hamil- 
tonian with quadratic terms A*<>Lc_ fc , has been widely 
accepted as a theoretical model for superconductors, 



where is a mean field of pair annihilation amplitude 
(c/j|C_fc|) and Cfc CT are annihilation operators of Fermion. 
The mean field violates the U(l) gauge symmetry, i.e., 
the total particle number is not preserved but fluctuat- 
ing. This charge fluctuation turns out to diverge in the 
thermodynamic limit @]. 

The most simple candidate of U(l) symmetry breaking 
terms is A&c|. + h.c, which implies injection of electrons 
and holes from out of the system. The concept of this 
term is directly connected to local charge fluctuation or 
doping. Although the previous study is limited to free 
fermion, such a ID system defined as H = J2k kc\ck + 
^Sfe( c I + c fe) nas been studied and was solved with a 
Jordan- Wigner transformation);! and with a canonical 
transformation^^- To clarify the meaning of A, let us 
describe the procedure of Ref.Q in detail; this spinless 
fermion chain can be mapped to semi-infinite XY model 
with a local magnetic field at the boundary. This local 
magnetic field in the xy plane turns out to be A. In Rcfs. 
0|T?J. A terms were introduced in different contexts. It 
is a common concept that the Hamiltonian is an effective 
one after tracing out of some environment. 

There are two motivations of the present work. The 
first is to clarify the properties of strongly correlated sys- 
tems with nonuniform charge fluctuation. Especially, a 
particle-hole symmetric Hubbard model with A will make 
it possible to realize the ground state which is a lin- 
ear combination of half-filled Mott insulating state with 
electron-doped and hole-doped states, which evokes the 
gossamer theory. It is interesting to evaluate fluctua- 
tion of total particle number as a direct effect of the 
U(l) symmetry breaking. The other is to test a mean- 
field type approach for the interchain hopping of fermion 
chains. When we consider a decoupling of interchain hop- 
ping t±cl a c±i a into Aio-cjg. with the mean-field type ap- 
proximation, the effective fields A^ can be identified as 
(t i cvr).|22j| In this approach, the Hamiltonian is consid- 
ered as an effective one obtained after tracing out neigh- 
boring chains in the quasi-one-dimensional (quasi-lD) 
systems. 

Since the Hubbard model has Coulomb interaction, we 
adopt the Density Matrix Renormalized Group (DMRG) 
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method, which is one of p owerful numerical methods for 
ID quantum systems [Hi HI HI- The DMRG method 
has achieved high precision in various ID systems while 
application to two-dimensional lattice is considered to be 
difficult. As an application to higher dimensions, we note 
that the DMRG method has already applied to quasi-lD 
spin systems with interchain couplings as mean fields in 
the manner of Ref. 14]. This paper may be first attempt 
to quasi-lD fermionic systems with interchain hopping 
as mean fields. 

The paper is organized as follows. In §[H]we construct 
a Hamiltonian with a generalized U(l) breaking term 
and mention a "bath" site introduced by the canonical 
transformation[I?ij . In 31111 we describe an application 
of the DMRG method to the Hamiltonian which does 
not conserve particle number. In iflVI we demonstrate 
the mean-field type approach for quasi-lD strongly cor- 
related electron systems. Finally we conclude the results 
and discuss problems of the mean-field type approach. 
In Appendix, the note for the canonical transformation 
is summarized. 



II. CORRELATED ELECTRON SYSTEMS 
WITH U(l) SYMMETRY BREAKING TERM 

Let us define a Hamiltonian with the generalized sym- 
metry breaking term term Ha 

H = H Q + H A , (1) 

where Hq can be one of the Hamiltonians for correlated 
electrons systems. In this paper we will restrict Hq to 
the Hubbard model defined as 

L-l 
i—1 (7 




where ct a are fermion operators and L is the system size. 
The symmetry breaking term Ha is defined as 

H A = J2 A trCi<r + A ia cl (3) 

where Ai a are considered as external fields at this stage. 
This model with nonzero Ai a breaks the particle num- 
ber conservation because Ha is not commutable with 
the total number of particles: [HA,N tot ] ^ 0, where 
Ntot = ^2cl a Ci a . This comes from the fact that Ha 
breaks the U(l) symmetry, where the global U(l) rota- 
tion defined as Ci a — > e l9 Ci a . 

Generally, external fields break some symmetry as 
magnetic fields break spin rotational symmetry. In ad- 
dition to the t/(l) symmetry-breaking, Ha also breaks 
the SU (2) spin-rotational symmetry, while the Hubbard 
model without magnetic field is an SU(2) invariant. The 



global SU(2) rotation is defined as Ci — > f/cj, where 
t Ci =* (cif,Ci|) and det U = 1. The symmetry break- 
ing term transforms under the SU (2) rotation as 

ff A (Ai) = ^c]A i + A| Ci 

i 

-^H' A = H A (U^A t ), 

where 'A, = (Aj-|-,Aj|). The symmetry is recovered 
only if Ai a — for all i, a. 



A. Particle-Hole Symmetry 

Let us suppose that Hq satisfies the particle-hole sym- 
metry, i.e., Hq is invariant under the usual particle- 
hole transformation on the tight-binding model: Ci a — > 
(— l) % (y icr . In other words, Hq is commute with a anti- 
unitary operator gHH, defined as 6 = KU v h, where K 
is a complex conjugation and C/ p h is the unitary operator 

defined as U p h = ^ L ]J ta (cia + (~1)' C L)- ^ satisfies 

e-^e = u-i Cia u ph = (-lfcl 

and one can show 

e-^e = Hq, 

e-^A^e = ^((-r/A,-), 

where t can be complex here and will be fixed t = 1 as the 
unit of energy in numerical calculations. The symmetry- 
breaking term Ha with A^ = (-l)'Ai is also invariant 
under the particle-hole transformation. Then, one can 
construct a Hamiltonian H = Hq + H a which preserves 
the particle-hole symmetry but breaks U(l) and SU(2) 
symmetries. It is easy to show total number of electrons 
of this Hamiltonian is half filled when the ground state is 
unique. The proof is following: because of the particle- 
hole symmetry _1 i/G = H, 9 _1 |gs) is also the ground 
state: H<d~ 1 \gs) = EggO^ 1 \gs). Since the ground state 
is unique, \gs) is equal to 9~ 1 |<?s) except for a phase 
factor. Then one gets 

(gs\N tot \gs) = (^le^A^ei^}, (4) 

and TVtot satisfies 

e^A^e = 2L-N tot . (5) 

From eqs.lQJ and JSJ, it is deduced that total number of 
electrons is half filled, (N to t) = L. It might be interest- 
ing to remind the reader that the half-filled Hamiltonian 
with the U(l) symmetry breaking term has some anal- 
ogy to a finite double occupancy in the half-filled case 
discussed in the gossamer superconducting theory pro- 
posed by LaughlinQ. 



B. Hidden even-odd Parity conservation 

To handle the fermion sign especially by the DMRG 
which will be discussed in the next section, let us con- 
sider a following extension of the Hilbert space by the 
canonical transformation |10| : 

where a is additional annihilation operator of a spinless 
fermion and satisfies 

{a, c icr } = 0, {a\c ia } = Q. 

Canonicality of this transformation can be easily shown 
as {ci a ,Ci'a'} — ~{cia,Ci' a '} = and {cio-jcj^,} = 

{cia, c i'cr'} ~ i^m'^ctct' • Moreover, when we denote the 
fermion operators in Ho as Ho({ci a }), Hq is invariant 
under this canonical transformation: 

H Q = H Q ({c ia }) -> So = -ffo({c iff }) = H Q ({c ia }) 

because Ho is made of invariant operators like (y ia Ci> a i — 
c \o- c i'cr'' The term H&, however, is modified as below, 

H A -» ^ = ^A* (T (a + a t )c l(T + A i(T cJ (T (a + a t ). 



This formula implies that a site is a "environment 
bath" site in the spirit of the Dynamical Mean Field 
Theory (DMFT) [Til. Since single fermion operators in 
H become bilinear, H conserves a even-odd parity of the 
total particle number defined below. 

The operator of the total particle number N tot is in- 
variant under the transformation, i.e., N to t = N to t- The 
total Hamiltonian H does not conserve -/V to t J}or the to- 
tal particle number including alpha N a = iVtot + a - 
That is [H, N a ] ^ 0. We note that such operators with 
a as N a has no correspondent operator, for example, N a 
which was defined before the transformation. 

The parity operator of N a is defined as 



(a) 



O) I"- 



(6) 



and satisfies P*P =J1 and P = . One can show the 
even-odd parity of N a is conserved, i.e., 

[H,P}=0. 

Since P is conserved, one can take simultaneous eigen 
states of H and P as 

H\*(E,p)) = E\*(E,p)), 
P\9(E,p)) = p\9{E,p)}, 

where E is eigen energy and p is ±1, because P 2 is 
identity as a operator and P is unitary. The ground 



IL„> IClr> IC^> IR„> 



(b) 




IL„> IClr> IC*> 



FIG. 1: The system in first iterative procedure for infinite sys- 
tem algorithm of the DMRG with the bath site "a" . The sys- 
tem size is enlarged from (a) L — 4 to (b) L = 6 . The a site, 
left block, center-left site, center-right site and right block are 
represented by indices \a n ), \L n ), \C„), \C„) and | -Re- 



states \^(E,p)) for H is doubly degenerated if the ground 
state \ty(E)) for H is unique. As shown in Appendix. 
1X1 one can show that expectation value of an arbi- 
trary operator O can be written as (E)\0\^ (E)) — 
($(E;+)\d\y{E;+)) = -)\6\*{E; -)). This 

means that any expectation value for the system H can 
be obtained in the system H . 



III. METHOD 

To study the present ID strongly correlated system 
without conservation of the total particle number, the 
DMRG method is used. As we implied in the previous 
section, the conservation of even-odd parity of particle 
number is required to calculate the fermion sign in the 
DMRG algorithm. In this section we will illustrate detail 
of the DMRG algorithm and mention the fermion sign. 
We note that the DMRG method has been applied to the 
different model which does not conserve total number of 
particles but conserves it's parity, called the BCS pairing 
Hamiltonian [l7| . 

First of all, we describe the iterative procedure of the 
DMRG for the Hamiltonian H with a site. Figure 
illustrates systems in first iterative procedure with en- 
larging the system size from L = 4 to L = 6. As seen 
in fig. ^b), the hopping terms between a site and each 
sites in H& become long range in the successive elonga- 
tion. Generally speaking, long-range hopping terms such 
as c\cl enlarge a numerical error but the numerical error 
is reduced because the a site is not renormalized in the 
iterative procedure. 

Next, to calculate the fermionic system, one should 
take care of the fermion signal- When the two lo- 
cal operators A, B are represented in bases \A n ),\B n ), 
one can usually get it's product (A n B m \AB\A n ' B m >) — 



4 



±(A n \A\A n i)(B m \B\B m i) , where the coefficient ± is the 
fermion sign. This formula is valid if the states, \A n ) and 
\B n ), have a fixed even-odd parity of particle number and 
operators A, B conserve the even-odd parity. Otherwise, 
states or operators are modified. For example such states 
as (l+c^jO) may change to (1 — c+)|0) when the states get 
the fermion sign. This change is hard to followed in the 
DMRG procedure. This is the reason why conservation 
of even-odd parity is required by the DMRG. 

The canonical transformation makes it possible to cal- 
culate the fermion sign because all operators conserve the 
even-odd parity as shown in ea. (|Al|) in Appendix. lAl 
Moreover, local bases, \L n ), \C„), \C„), \R n ) and \a n ), 
have fixed parity as illustrated below. 

Then, all local bases in each steps should have fixed 
even-odd parity of particle number. We explain this with 
one step of the iterative procedure below. 

1. Here we suppose local bases have fixed parity of 
particle number as in L = 4 system. 

2. Make the matrix elements of the Hamiltonian H 
from local operators represented by local bases, 
taking care of the fermion sign. The Hamiltonian 
H is block diagonalized and has even and odd par- 
ity sector. 

3. Calculate the ground states \^(E gs ,p)) and it's en- 
ergy E gs , where p — ± is even-odd parity of particle 
number. 

4. Make the matrix elements of four density matrices 
Pl/r,±; 

Pl,± = Tr p±, p R ± = Tr p±, (7) 

R,C R ,a L,C L ,a 

where p± = \®{E gs , ±))(* (E gs , ±)|. One can 
show that these density matrices are block di- 
agonalized into even and odd parity sector, i.e., 
[Pl/r,±, e 17rJVl /«] = 0, as the Hamiltonian is. We 
note that if the Hamiltonian conserves particle 
number, [H, N t ot] = 0, the density matrices for 
the left and right block conserve particle number, 
[Pl/r, N L/R ] = 0. 

5. Diagonalize Pl/r,± and select the lowest m eigen- 
values and their eigenvectors called renormalized 
bases, which have fixed parity of particle number 
because Pl/r,± conserve the parity. 

6. Remake matrix elements of all local operators in 
the renormalized bases. Then, renormalized bases 
are the next local bases and satisfy the supposition 
in the step 1. 

As described in the procedure, since H conserves even- 
odd parity, one can show that each local base in each 
steps has fixed parity of particle number. That is, one 
can calculate the fermion sign. 



Finally, we note that the number of the states of left 
and right block is used up to about 60 and the truncation 
error is less than 10 -4 in the following results. Since we 
deal with general Ai a which depends on the site, the 
DMRG method for the random system^} is employed. 

The expectation value is deduced as 

( A ) = lj2^ E ^P)\AmE gs ,p)), (8) 
p=± 

in order to avoid the numerical error. 



IV. RESULTS 

In the following, we will show results of a total-charge 
fluctuation and a particle number in ijlV Al which are the 
direct effect of existence of external fields Aj CT . In 3IVBI 
we will show tentative demonstration of the mean-field 
theory to deal with the interchain hopping. 

We note that we restrict the system size L to even 
to obtain the unique ground state \^f(E)) for H. In 
the DMRG calculation, H is used and the ground states 
\^/(E,p)) for H is doubly degenerated. This degeneracy 
is artificial due to the canonical transformation. 

A. Total-charge fluctuation 

When the Hamiltonian Hq is commutable with N, the 
ground state for Hq has a fixed number of electrons and 

the total-charge fluctuation, AN = J (N 2 ) - (N) 2 , is 
zero when the ground state is unique. On the other 
hand, since external fields breaks the U (1) symmetry, 
the total-charge fluctuation becomes finite. The ground 
state for nonzero Ai a becomes a linear combination of 
electron-doped states and hole-doped states. It implies, 
even at half-filling, "effective carriers " are introduced by 
the nonzero A^. 

Since large charge gap prefers no fluctuation, when the 
Coulomb interaction U becomes infinite, doped states 
are not allowed at half-filling and AN becomes zero. In 
FigiU total charge fluctuation AA" as a function of 1/U 
is plotted as an example of the particle-hole symmetric 
Ai a . Infinite U gives no fluctuation AN — and finite 
total-charge fluctuation is proportional to 1/U. This fig- 
ure is only for the small system and relatively large Ai a 
but the qualitative character that AN is proportional to 
1/U is general. 

It should be noted that maximum of total-charge fluc- 
tuation AN is order 1 as fig[21 shows AA~ ~ 0(1) in 
large 1/U region. That is, the present model can not 
reproduce the BCS ground state with AN ~ 0{L x l 2 ). 
This is important for the charge compressibility defined 
as k(p) = lp = j- ^gisil _ Let us summarize this prop- 
erty for the Mott insulator and the metal: in the Mott 
insulating state, the charge compressibility diverges at 
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FIG. 2: An example of total charge fluctuation AN as a 
function of 1/U for system-size L — 4 at half- filling with 
the Hamiltonian H = Ho + Ha at t = 1 and particle-hole 
symmetric A nCT = cos(n7r/4), where {N to t) = L is satisfied 
numerically. This figure is only for the small system and rel- 
atively large Aj CT but the qualitative character that AiV is 
proportional to 1/U is general. 



fi = ±A c /2 and K(/i]_= at — A c /2 < n < A c /2, where 
A c is the charge gapl9j. In the metallic case, the charge 
compressibility is proportional to the density of states at 
the Fermi energy. In non-zero Ai a case, the charge com- 
pressibility at half-filling becomes finite but is expected 
to become zero at fi = in the limit L — > oo because the 
magnitude of A7V is order 0(1). 



B. inter chain hopping as mean fields 

In this subsection, let us consider the interchain hop- 
ping of ID Mott insulators. When we take an ansatz 
of the mean field type, external fields A,:„ can be de- 



termined self-consistently as — t±c 
adopt the self-consistent equation: 



We 



FIG. 3: Converged Ai a = Aj as a real function of site i with 
varying U = 0,3,4. The system size is L = 20. 




FIG. 4: Stabilization energy AE = E(0) - E(A irT ) as a func- 
tion of t± with converged Ai a for U/t — 2. Points are ob- 
tained by the DMRG and lines are fitted with a linear func- 
tion. 



A ia = -t±(c icr ). 

Although the interchain hopping also gives rise to the 
effect of the band structure and the dimensionality, the 
effect of the band structure is not taken into account in 
this approach. The meaning of tj_ in the self-consistent 
equation is the strength of charge fluctuation in the per- 
pendicular direction with the general band structure. 

In the DMRG method we used the transformed Hamil- 
tonian, H, and the transformed self-consistent equation 
as A icr = — t±((a + a' We note that A ia is lim- 
ited to real number for simplicity. In Fig|3] some result 
of converged A^ after the self-consistent loop for L = 20 
are plotted. Ai a converged smaller value as U increased, 
which implies a quantum phase transition from nonzero 
Ai„ to zero A icr . 



To clarify the transition, one can define the stabiliza- 
tion energy AE = E(0) — E(Ai a ) 1 which means the en- 
ergy gain due to charge fluctuation in the perpendicular 
direction, where E(Ai a ) is the ground state energy with 
converged Aj CT . There are two simple limits: infinite U 
limit and small t± limit. In both case A^ converged to 
about zero and the ID Mott insulator is realized. In Fig0] 
the stabilization energy is plotted as a function of tj_ . In- 
creasing ij_ means that length between chains is changing 
more closely, which corresponds to applying pressure. In 
small tj_ region, AE becomes zero, which is identified 
as the ID Mott insulator phase. Actually, converged A ia 
and AiV" are zero there. There is a transition from the ID 
Mott insulator phase to the symmetry breaking phase as 
t± increase. Extrapolated points in FigQJshow that the 
transition point t\/t is about 0.16 for U/t — 2. Since the 
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charge gap of Mott chain of U/t = 2 is A c = O.170, the 
naive criterion ij_ ~ A„|20| is reasonable in this analysis. 



V. CONCLUSION 

In conclusion, we have studied the symmetry breaking 
of U (1) charge and SU (2) spin due to external fields A ia , 
which lead to nonzero total-charge fluctuation AN. Fi- 
nite AN means the coherent ground state is a linear com- 
bination of electron-doped states and hole doped states. 
We have applied the DMRG method to the particle-hole 
symmetric Hubbard chain with Ai a and have actually 
demonstrated that the total-charge fluctuation AN at 
zero-temperature is linear with respect to l/U even at 
half- filled case (N to t) = L. 

Considering Ai a as a mean field of interchain-hopping 
tentatively, we have obtained the quantum phase transi- 
tion from the ID Mott insulator to the symmetry break- 
ing phase as t± increases, i.e., pressure increases. The 
transition point t c ± in Figgis close to the charge gap A c . 
In the symmetry breaking phase, effectively doped car- 
riers itinerate between chains because of non-zero A^. 
Since the difference between two phases is whether the 
interchain hopping becomes relevant or not, one may say 
this transition deconfinement transition|20j. However, 
since this approach is limited to small A.^, clarification 
of the symmetry breaking phase is remained as future 
work. The clarification is interesting work because it is 
well known that destroying the Mott phase by apply- 
ing the pressure, i.e., increasing t±, is typical for high-T c 
cuprates. 

As described in flIVI the magnitude of AN as a func- 
tion of the system size is constant while the BCS theory 
gives AN ~ (DfX 1 / 2 ). This property may be related to 
the fact that we dropped the anti-commutation relation 
between Aj CT and fermion operators in the Hamiltonian. 
That is, expectation value (cj CT ) was a fermionic operator 
before taking the average as a mean-field. As pointed out 
in ref. |2l|, this fact gives the limitation of this "mean 
field" approach. We can find the another way to bosonize 
the mean field and this result will be reported in future 
work. 

As a technical outlook, we have used only the infinite 
method of the DMRG in the self-consistent calculation. 
The combination of finite method and self-consistent loop 
will improve the cost of calculation time, where self- 
consistent field is calculated at center block in sweep of 
the finite method. In this method, DMRG is combined 
with the mean-field theory more closely. 
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APPENDIX A: CANONICAL 
TRANSFORMATION 

To summarize properties of the canonical transforma- 
tion, we introduce the majorana fcrmions defined as 

a + — a + 

a- = — i(q — a^), 

which are unitary and Hermite and satisfy anti- 
commutation relations {a+,a_} = 0, {a±,Ci a } = 0. It 
can be proved that even-parity of P defined in eq. © is 
anti-commutable with a_: {P, a_} = 0. In addition, 
the canonical transformation changes any operator for 
the system H, denoted by O, into O which satisfies 

[P,d] = [a-,d] = Q, (Al) 

because O consists of c,- CT , c[ a and a+, not of a_ and the 
operator as a polynomial function of fermion operators 
does not have terms of odd degree. Finally, since the 
transformation is canonical and the new vacuum is de- 
fined as Ci CT |0) = a\0) = 0, it can be shown that 

(0|O|0) = (0\O\0). (A2) 

To define the canonical transformation of the states, 
we define bases explicitly as 

2L 

\r>~\{m}) = n^) ni l°>- ( A3 ) 

i=l 

After the canonical transformation, the bases are written 
as 

2L 

\f)--=\{m}) = H(a + 4r\o). ( A4 ) 

i=i 

Since a\ = 1, the bases has even-parity of P: P\I) = \I), 
where P is defined in eq.lJBJ. When we define |J;+) = 
| J) and 1 7; — ) = a_|7), one can easily show that a set 
of 2 x 4 L bases \I; ±) is the ortho-normalized complete 
set and the bases satisfy P|7;±) = ±|7;±) because of 
{P,a-} = 0. 

With ca. (|Al|) and a 2 _ = 1, one can show 
(T;+\d\T';+) = (T;-\d\T;-). With eq.^ one can 
also show (7|C|7) = (7; +\0\T\ +). Combining them, we 
summarize 

(7|0|7') = (7;+|5|7';+) = (7;-|a|7';-). (A5) 

This means the block diagonalized operator has the same 
matrix element for even-odd sectors. Since equation (|A5|) 
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is true when O is the Hamiltonian H, the eigen vectors 
of H and H can be the same: 

|*(£)> = £Cr(£0|i) (A6) 

7 

\*(E,±)) = £ 0(25)1/; ±), (A7) 
where elements Ci(E) satisfy (I\H\I')Cr{E) = ECi(E). 



We note that degenerated eigen vectors \^(E, ±)) satisfy 
($(E, P )\y(E',p'))=6 E E>S PP >. 

We conclude that the expectation value for any oper- 
ator O can be written as 

(*(£)|O|*(£0) = <$(25;+)|<5|$(25;+)) 

= <$(25;-)|0|$(25;-)). (A8) 
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